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ABSTRACT 

The first stars form in dark matter halos of masses ^IQ^Mq as suggested by an increasing number of 
numerical simulations. Radiation feedback from these stars expels most of the gas from their shallow 
potential well of their surrounding dark matter halos. We use cosmological adaptive mesh refinement 
simulations that include self-consistent Population III star formation and feedback to examine the 
properties of assembling early dwarf galaxies. Accurate radiative transport is modelled with adaptive 
ray tracing. We include supernova explosions and follow the metal enrichment of the intergalactic 
medium. The calculations focus on the formation of several dwarf galaxies and their progenitors. In 
these halos, baryon fractions in 10® Mq halos decrease by a factor of 2 with stellar feedback and by 
a factor of 3 with supernova explosions. We find that radiation feedback and supernova explosions 
increase gaseous spin parameters up to a factor of 4 and vary with time. Stellar feedback, supernova 
explosions, and H2 cooling create a complex, multi-phase interstellar medium whose densities and 
temperatures can span up to 6 orders of magnitude at a given radius. The pair-instability supernovae 
of Population III stars alone enrich the halos with virial temperatures of 10^ K to approximately 
10~^ of solar metallicity. We find that 40% of the heavy elements resides in the intergalactic medium 
(IGM) at the end of our calculations. The highest metallicity gas exists in supernova remnants and 
very dilute regions of the IGM. 

Subject headings: cosmology: theory — galaxies: dwarf — galaxies: high-redshift — stars: formation 



1. MOTIVATION 

The majority of galaxies in the universe are low- 
luminosity, have masses of ^10® solar mas s es, and are 
known as d warf galaxies (jSchechteil 119761: lEllisI 119971 : 
lMateol[l998 h. Galaxies form hierarchically through nu- 
mero u s mergers of smaller c ounterparts ( (Peebles fc Dicke 
119681 IWhite fc ReesI Il978[ ). whose properties will in- 
evitably influence their parent galaxy. Dwarf galaxies 
are the smallest galactic building blocks, and this leads 
to the question on even smaller scales: how were dwarf 
galaxies influenced by their progenitors? 

A subset of dwarf galaxies, dwarf sphcroidals 
(dSph), have the high e st ma s s-to-l ight ratios 
(|de Blok & McGaue-hl [l997t iMateol [l998h and con- 
tain a population of meta l-poor sta rs that arc similar 
to G alactic halo stars (Tol stov et al . 2004; Hclm i~et al.l 
[20061) . There is a metallicity floor exists of 10 ^ and 
10~^ of so lar metallicity in dSph and halo stars, re - 
spectively (|Beers fc Christliebl l2005t iHelmi et aLll2006D . 
Stellar metallicities increase with time as previous 
stars continually enrich the interstellar medium (ISM). 
Hence the lowest metallicity stars are some of the oldest 
stars in the system and can shed light on the initial 
formation of dwarf galaxies. This metallicity floor also 
suggests that metal enrichment was widespread in dark 
matt er halos befor e low-m ass stars could have formed 
(e.g. iRicotti et all l2002b| ). Supernovae (SNe) from 
metal-free (Pop III) stars generate the first metals in 
the universe and may supply the necessary me tallicity 
to form the most metal-poor stars observed (jFerraral 
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[l998HMadau et al.l[200H iNorman et al.ll2004f ). 

Dwarf galaxy formation can be further constrained 
with observations that probe reionization and semi- 
analytic models. Observations of luminous quasars 
powered b y supermassive blac k holes (SMBH) o f mass 
^O^Mq (iBecker et al.l [200l iFan et all l200l. [200il 



and low-luminosity galaxies fHu et al. '"'2002'; 'lye et ajj 
2006; .Kashikawa et aL 2006; Bouwcns & lUingwort^ 
2006t IStark et all l2007l ) at and above redshift 6 indi- 



cate that active star and BH formation began long 
before this epoch. Semi-analytic models have ar- 
gued that cosmological reionization was larg ely caused 
by l ow-luminosity dwarf gala xies ( JIaiman 



Ce3[2003i: SomcrviUc & Livioi 1200.1 IWise & Abe! 



1997 



2005 



Haiman fc Brvan .2006). Some of the most relevant pa- 



rameters in these models control star formation rates, 
ionizing photon escape fractions, metal enrichment, and 
the minimum mass of a star forming halos. They are 
usually constrained using (i) the cosmic microwave back- 
ground (CMB) polarization observation from Wilkin- 
son Microwave Anisotropy Probe (WMAP) that mea- 
sures the optical depth of electron scattering to the 
CMB (P age et al.ll2007[ ). (ii) Gunn-Peterson troughs in 
z ~ 6 quasars, and (iii) numerical simulations that 
examine neg ative and positiye feedback of radiation 
backg rounds (iMachacek et al.l l2001l. 120031: lYoshida et all 



120031 iMesigner et al.l |2006[) . Radiation hydrodynami- 
cal ab initio simulatio ns of the first stars (lYoshida et al.l 
l2007HAbel et al.ll2007f ) and galaxies can further constrain 
the parameters used in semi-analytic models by analyz- 
ing the impact of stellar feedback on star formation rates 
and the propagation of H II regions in the early universe. 
Moreover, these simulations contain a wealth of informa- 
tion pertaining to the properties of Pop III star forming 
halos and early dwarf galaxies that can increase our un- 
derstanding of the first stages of galaxy formation. 
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First we need to consider Pop III stars, which form 
in the progenitor halos of the first galaxies, to cap- 
ture the initial properties of dwarf galaxies. Cosmo- 
logical numerical studies have shown that massive (30- 
300 Mq) Pop III stars form in dark matter halos with 
masses ^lO^M^ (lAbel et al.ll2003 iBromm et all 120021: 



Yoshida et all2006bl:lGao et al.ll2007i : lO'Shea fc Normanl 
200 /t ). Recently. lYoshida et all and iTurk et all (|2008D 



followed the gaseous collapse of a molecular cloud that 
will host a Pop III star to cosmologically high number 
densities of 10"'^^ and 10^^ cm~^, respectively. The for- 
mer group thoroughly analyzed the gas dynamics, cool- 
ing, and stability of this free-fall collapse. The lat- 
ter group observed a protostellar core forming with 10 
Jupiter masses that is bounded by a highly asymmetric 
protostellar shock. Both groups found no fragmentation 
in the fully molecular core that collap ses into a single, 
massiv e '^IOOMq star. Furthermore, lOmukai fc Pallal 
(|2003! ) determined that accretion may halt at the same 
mass scale, using protostellar models even for different 
mass accretion histories. 

Pop III stars with stellar masses roughly between 
140 and 260 M© end their life in a pair-instability SN 
(PISN) that releases 10^^ — 10^'^ ergs of energy and 
tens of solar masses of heavy elements into the ambi- 
ent medium (Barkat, Rakayy , fc Sack 1967; Bon d et al.l 
119841: iHeger fc Wooslevl l2002[ ). These explosions are an 
order of magnit ude larger than typical Ty pe II SNe in 
both quantities (jWooslev fc Weaverl Il986f ). such explo- 
sions energies are larger than the binding energies in 
their low-mass hosts, e.g., 2.8 x 10^° ergs for a lO^M© 
halo at redshift 20. Thus gas structures in the host halo 
are totally disrupted and expelled, effectively enriching 
the surrou nding intergalactic me d ium (IGM) with the 
SN e jecta (B romm fc Loebl [200I iKitavama fc Yoshidal 
120051 [Grcif et al. 2007). The combination of the shal- 
low potential well and large explosion energy suggests 
that these events are good candidates for enriching the 
first galaxies and IGM. Outside of the pair-instability 
mass ran ge. Pop II I stars die by directly collapsing into 
a BH ( Hc ger et al.l [2003). possibly providing the seeds of 
high-redshift quasars in galaxies th at are associated with 
the rarest density fluctuations (e. g . ^Mad au fc Reed 120011 : 
IVolonteri et 1111120051 iTrenti fc St iavelh 2 00l. 

One-dim ensional calcula tions ([W halcn et _al.l l2004t 
iKitavarnaT ct al. 2004; Kita vama fc Yo shida 20Q1) and re- 
cent three-di mensional rad i ation hydrqdyna mical simula- 
tions (Yoshida et al.ll2007t lAbel et all 12007 ) have inves- 
tigated how the Pop III stellar feedback affects its host 
halo and nearby cosmic structure. In addition to SNe, 
H II regions surroun din g Pop III stars, which have lumi- 
nosities --IO^Lq (|Schaereiil2002i) . alone can dynamically 
affect gas at distances up to a fe w proper kp c. Ioniza- 
tion fronts and H II regions fsee fY^TO [19861 for a re- 
view) have been extensively studied in literature on star 
formation since Stromgren (1939). Stellar radiation gen- 
erates an ionization front that begins as a R-type front 
and transforms into a D-type front when its speed slows 
to twice the sound speed of the ionized gas. Then a 
strong shock wave forms at the front and recedes from 
the star at ~30 km s~^. The ionization front decouples 
from the shock wave and creates a final H II region that 
is 1 - 3 proper kpc in radius for massive Pop III stars 
residing in low-mass halos. The ionized gas is warm 



{'^3 X 10'* K) and diffuse (~1 cm~^). The shock wave 
continues to accumulate gas and advance after the star 
dies. Eventually it stalls in the IGM, but in the pro- 
cess, it r educes the bary on fraction of the halo below one 
percent (jYoshida et allboOTtlAbel et al.l l2007^. 

Clearly the number of progenitors of a given galaxy as 
well as the star formation and feedback history of the 
progenitors will play a role in shaping all of its proper- 
ties. But how much? If most stars of a galaxy are formed 
later, will the earliest episodes not be entirely negligible? 
To start addressing these questions, we have carried out 
a suite of simulations that include accurate three dimen- 
sional radiative transfer and the SN explosions of Pop 
III stars and have followed the buildup of several dwarf 
galaxies from those Pop III star hosting progenitors. The 
Pop III radiative and SN feedback dramatically alter the 
properties of high redshift dwarf galaxies, and we dis- 
cuss some of the most striking differences here. We leave 
a more detailed exposition of star formation rates, star 
forming environments, and the beginning of cosmic reion- 
ization for a later paper. 

In the following section, we detail our cosmological, ra- 
diation hydrodynamics simulations and the star forma- 
tion algorithm. Then we describe the global characteris- 
tics of dwarf galaxies that forms in our simulations in |j3l 
There we also focus on metal enrichment of star forming 
halos and the IGM, arising from PISNe. In fjll we discuss 
the implications of our findings on the paradigm of high- 
redshift galaxy formation by including II2 chemistry and 
Pop III star formation and feedback. We summarize in 
the last section. 

2. RADIATION HYDRODYNAMICAL 
SIMULATIONS 

We use the Eulerian AMR hy drodynamic code 
Enzo (jBrvan fc Normanl Il997l Il999f i to study the im- 
portance of primordial stellar feedback in early galaxy 
format ion. Enzo uses a n n-body adaptive particle-mesh 
solver ()Couchmanlll991l ) to follow the dark matter (DM) 
dynamics. We first describe the setup of our simulations. 
We then detail our star formation recipe for primordial 
star formation. Also we have implemented adaptive ray 
tracing into Enzo whose description concludes this sec- 
tion. 

2.1. Simulation Setup 

We perform two cosmological realizations with differ- 
ent box sizes and random phases and WMAP 1 year pa- 
rameters of {h, r^ A, r^Mi o'R, n) = (0.72, 0.73, 0.27, 
0.024/1-2, Q 9^ ^-^ (iSoergel et al.lf200l . Table [1] summa- 
rizes the details of these simulations. The characteristics 
of the individual halos studied here are not affected by 
the signific antly different WM AP third year parameters 
(WMAP3: [Spergel et al.]|2007D . which do affect the sta- 
tistical properties of such halos. We have verified that 
nothing atypical occurs during the a ssembly of the halos 
studied here fsee lWise fc Abe]||2007l hereafter Paper I). 

The initial conditions are the same as in Paper I. 
To simplify the discussion, simulation A will always be 
quoted first with the value from simulation B in paren- 
theses. They both have a top grid with a resolution of 
128"^ with three nested subgrids with twice finer resolu- 
tion and are ini tialized at z = 129 (119) with the COS- 
MICS package ()Bertschingeijll995l . l200l[ ). The box size 
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TABLE 1 

Simulation Parameters 



Name 


I 

[Mpc] 


SF 


SNe 


Npart 




Ng,id 




N 


cell 


SimA-HHe 


1.0 


No 


No 


2.22 


X 


10^ 


40601 


1.20 


X 


10* (4943 ) 


SimA-RT 


1.0 


Yes 


No 


2.22 


X 


10^ 


44664 


1.19 


X 


10* (4933 ) 


SimB-HHe 


1.5 


No 


No 


1.26 


X 


107 


21409 


6.51 


X 


107 (4023) 


SimB-RT 


1.5 


Yes 


No 


1.26 


X 


10^ 


24013 


6.54 


X 


107 (403^) 


SimB-SNe 


1.5 


Yes 


Yes 


1.26 


X 


107 


24996 


6.39 


X 


107 (4003) 



Note. — Col. (1): Simulation name. Col. (2): Box size. Col. (3): Star 
formation. Col. (4): Supernova feedback. Col. (5): Number of dark matter 
particles. Col. (6): Number of AMR grids. Col. (7): Number of unique grid 
cells. 



is 1.0 (1.5) comoving Mpc. The innermost grid has an 
effective resolution of 1024^ with DM particle masses of 
30 (101) Moand a side length of 250 (300) comoving kpc. 

Regions of the simulation grid arc refined by two when 
one or more of the following conditions are met: (1) 
baryon density is greater than 3 times QbPoN^''^~^'^\ (2) 
DM density is greater than 3 times OcdmPqA^'^^^'^'': and 
(3) the local Jeans length is less than 16 cell widths. Here 
= 2 is the refinement factor; I is the AMR refine- 
ment level; 4> = —0.3 causes more frequent refinement 
with increasing AMR levels, i.e. super-Lagrangian be- 
havior; po = 3Hq/8ttG is the critical density; and the 

Jeans length, Lj = yj 15kT / Att pG iimn , where Hq, k, T, 
p, p,, and ruH are the Hubble constant, Boltzmann con- 
stant, temperature, gas density, mean molecular weight 
in units of the proton mass, and hydrogen mass, respec- 
tively. The Jeans length refinement insures that we meet 
the Truelove criterion, which requires the Jeans length to 
be re solved by at least 4 cells on each axis (| Truelove et al.l 
Il997f) . Further refinement is only allowed in the initial 
innermost grid that has a comoving side length of 250 
(300) kpc. We enforce a maximum AMR level of 12 in 
these simulation that corresponds to a resolution limit of 
2.9 (1.9) comoving parsecs. 

We use the nine species (H, H+, He, He+, He++, 
e~ , H 2, Hj", H~) non-equilibrium chemistry model in 
Enzo (|Abel et all Il997t lAnninos et all I1997D and the 
H2 cooling rates from iGalH fc Fallal (|1998D . Compton 
cooling and heating of free electrons by the CMB and 
radiative losses from atomic and molecular cooling are 
computed in the optically thin limit. 

We focus on the region containing the most massive 
halo in the simulation box. We perform three calcula- 
tions - simulation A with star formation and radiation 
transport (RT; SimA-RT), simulation B with star for- 
mation and RT (SimB-RT), and simulation B with star 
formation, RT, and SNe (SimB-SN). We end the calcula- 
tions at the same redshift the halo with a virial tempera- 
ture of 10^ K collapses at z = 15.9 (16.8) in the hydrogen 
and helium cooling runs (HHe) of Paper I. 

2.2. Star Formation Recipe 

Star formation is mo delled through a n ext ension 
(lAbel et all [200I of the iCen fc OstrikeJ (fl99l algo- 
rithm that automatically forms a star particle when a 
grid cell has 

1. an overdensity exceeding 5 x 10^ 

2. a converging velocity field (V • v < 0) 



3. rapidly cooHng gas (icooi < ^dyn) 

4. an H2 fraction greater than 5 x lO^"'. 

Then we remove half of the gas from the grid cells in a 
sphere that contains twice the stellar mass, which is a 
free p arameter. Once these criteria are met, Abel et al.l 
(|2002f ) showed that a Pop HI star forms within 10 Myr. 
For this reason, we do not impose th e Jeans instability 
requirement used in ICen fc Ostrikeil and do not follow 
the collapses to stellar scales. We allow star formation 
to occur in the Lagrangian volume of the surrounding 
region out to three virial radii from the most massive 
halo at z = 10 in the dark matter only runs as discussed 
in Paper I. This volume that has a side length of 195 
(225) comoving kpc at z = 30 and 145 (160) comoving 
kpc at the end of the calculation. 

Runs with only star formation model all Pop III stars 
with M^ = 100 Mq that live for 2.7 Myr and emit 1.23 x 
10^° hydrogen ionizing photons per second. After its 
death, the star particle is converted into an inert 100 
Mq tracer particle. The SNe runs use M* = 170 Mq that 
results in a lifetime of 2.3 Myr and 2.57 x 10^° ionizing 
photons per seco nd, in accorda nce with the no mass loss 
stellar models of ISchaereil (|2002) . 

When a 170 Mq star dies, it injects Esn = 2.7 x 10^^ 
erg of thermal energy and SO/go Mq, of metals, appropr i- 
ate for a PISN of a 170 Mq star (|Heger fc Wooslevll2002D , 
into a sphere with radius rsN = 1 pc centered on the 
star's position. The mass contained in the star particle 
and associated metal ejecta are evenly distributed in this 
sphere. The mass of the star particle is changed to zero, 
and we track its position in order to determine the num- 
ber of stars associated with each halo. We also evenly 
deposit the SN energy in the sphere, which changes the 
specific energy by 



Po + PSN 



(1) 



where po and eo arc the original gas density and spe- 
cific energy, respectively. Here psN — Mi,/Vs^ is the 
ejecta density; esN = -EsnZ-M^/Vsn is the ejecta specific 
energy; Vsn is the volume of a sphere with radius rsN- 
In order not to create unrealistically strong shocks at 
the blast wave and for numerical stability reasons, we 
smoothly transition from this energy bubble to the am- 
bient medium, using the function 



/(r) = A <^ 0.5 - 0.5 tanh 



B[ 1 



(2) 
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Here A is a normalization factor that ensures J f{r)dr — 
1, and B controls the rate of transition to the ambient 
medium, where the transition is steeper with increasing 
B. We use A = 1.28 and _B = 10 in our calculations. 

We continue to use the nine-species chemistry model 
and do not consider the additional cooling from metal 
lines and dust. We follow the hydrodynamic transport 
of metals from the SNe to the enrichment of the sur- 
rounding IGM and halos. 

2.3. Radiative Transfer 

For point sources of radiation, ray tracing is an ac- 
curate method to calculate and evolve radiation fields. 
However, millions of rays must be cast in order to ob- 
tain adequate ray sampling at large radii. We use adap- 
tive ray tracing (jAbel fc Wandeltl l2002f l to overcome 
this d ilemma associated with ray tracing (cf. lAbel et ahl 
|2007( ). We initially cast 768 rays, i.e. level three in 
HEALPix (Gorski ct al. 2005), from the radiation source. 
The photons contained in the initial rays are equal, and 
their sum is the stellar luminosity. Rays are split into 
4 child rays, whose angles are calculated with the next 
HEALPix level, if their associated solid angle is greater 
than 20% of the cell area (Ax)^. Photons are distributed 
evenly among the children. This occurs if the ray travels 
to a large distance from its source, or the ray encoun- 
ters a highly refined AMR grid, in which adaptive ray 
tracing accurately samples and retains the fine structure 
contained in high resolution regions. 

The rays cast in these simulations have an energy Epi^ 
that is the mean energy of hydrogen ionizing photons 
from the stellar source. For 100 Mq and 170 Mq, this 
energy is roughly equal at 28.4 and 29.2 eV, respectively, 
due to the weak dependence of the surface temperature 
of primordial stars on stellar mass. The rays are trans- 
ported at the speed of light on timesteps equal to the hy- 
drodynamical timestep on the current finest AMR level. 

The radiation transport is coupled with the hydrody- 
namical, chemistry, and energy solvers of Enzo. Here we 
only consider hydrogen photoionization. We first calcu- 
late the photoionization and heating rates caused by each 
ray and then sub-cycle the chemistry and heating solvers 
with these additional rates on every radiation timestep. 
Next we advance the hydrodynamics of the system with 
the usual adaptive timesteps. 

The hydrogen photoionization rate is computed by 



^ riHi Vceii Aiph ' 

where Pq is the incoming number of photons, tihi is the 
number density of neutral hydrogen, Vccii is the volume 
of the computational grid cell, and r = nni ctri dl is the 
optical depth. Here am is the cross section of hydro- 
gen, and dl is the distance travelled by the ray through 
the cell. The heating rate is computed from the excess 
photon energies of the photoionizations by 

Pi = kph{Eph — Ei) 

= ^^iEp,-E.). (4) 

In the case of hydrogen ionizing photons, Ei = 13.6 eV. 
In each radiation timestep, the number of photons ab- 
sorbed, i.e. Po(l ~ e""^), is subtracted from the ray. The 



ray is eliminated once most of the associated photons 
(e.g. 99%) are absorbed or the ray encounters a highly 
optically thick region (e.g. r > 20). 

To model the H2 dissociating (Lyman- Werner; LW) 
radiation between 11.2 and 13.6 cV, we use an opti- 
cally thin 1 /r^ radia t ion fi eld with luminosities calcu- 
lated from ISchaereij (|2002D . We use the H2 photo- 
dissociation rate coefficient for the Solomon process of 
kdiss — 1-1 X lO^-FLw s~\ whe re Ft.w is the LW flux 
in units of erg s~^ cm"^ Hz"^ ()Abel et al.lll997[ ). We 
do not consider the self-shielding of LW photons. Be- 
cause molecular clouds only become optically thick in 
the late stages of collapse and above c olumn densities 
of '^10^^ cm~^ (jPraine fc Bertoldilfl996l) . we expect our 
results to not be drastically affected by neglecting LW 
self-shielding. Additionally, LW self-shielding may be 
unimportant up to column densities of 10^" — 10^^ cm~^ 
if the mediu m contains very large v elocity gradients and 
anisotropics (jGlover fc Brandll2001l ). 

In principle, the propagation of LW photons can be 
calculated with our ray tracing software. However, the 
computational expense of this additional coupled calcula- 
tion is large because the LW radiation field is essentially 
optically thin and must be traced through all grid cells 
in the simulation. In comparison, the hydrogen ioniz- 
ing photons are only traced within the H II region whose 
radii are at most 5 proper kpc. To test the validity of our 
approximation, we calculate the H2 optical depth from a 
currently shining star to a spherical surface with a radius 
of 3 proper kpc, which is plotted in Figure [1] The top 
panel depicts the effect of sclf-shiclding above a column 
density of lO^"* cm~^ and can be fitted with the power 

law, F = FqN~1''^. Here F and Fq are the transmitted 
and unobscured fiux, respectivel y, and A^i4 is the H2 opt i- 
cal depth in units of 10 cm" -2 ( Draine fc Bertoldill996D . 

We calculate the H2 optical depth in static data out- 
puts in SimA-RT, SimB-RT, and SimB-SNe with adap- 
tive ray tracing. The H2 abundances are taken from 
the simulations, which use the optically-thin assump- 
tion. This causes the optical depths in Figure [T] and 
self-shielding to be underestimated. Nonetheless, it is 
useful to estimate the amount of H2 self-shielding in 
these high redshift halos. The stars are located in dark 
matter halos with masses 3 x 10'', 2 x 10^, 10^, and 
5 X IQS Mq at z « (16.8, 19, 22, 30), respectively. In 
halos with M < IO^Mq, almost all of the molecular hy- 
drogen within the host and neighboring halos are not 
self-shielding, which starts to occur above column den- 
sities of lO^'' cm~^. Above this halo mass, a fraction 
of the H2 column densities rise above this critical value 
for H2 self-shielding. Thus at lower redshift, we could 
be underestimating the star formation rate as overdense 
clumps may be shielded to the LW radiation, especially 
in larger halos or when a primordial star is living in a 
neighboring massive halo. 

3. RESULTS 

In this section, we first discuss star formation in dwarf 
galaxy progenitors. Then we focus on the global charac- 
teristics of the most massive halo. We detail the different 
ISM phases. Metal transport from PISNe and the associ- 
ated metal-enriched star formation history are discussed 
last. 
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TABLE 2 

Global halo properties with star formation 
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log H2 Column Density [cm"^] 

Fig. 1. — Top: Normalized, trans mitte d Lyman- Werner flux as a 
function of H2 column density [ Draine & : Bcrtoldi 1996). Bottom: 
Histograms of H2 column density at 3 proper kpc from a currently 
shining star. The various lines are evaulated at different redshifts 
and host halo masses, which are approximately (30, 10^''^ Mq) for 
solid black, (22, W^^-^Mq) for dotted red, (19, 1O"^-3M0) for dashed 
green, and (16.8, IO^'^Mq) for dash-dotted blue. The vertical thin 
line at 10^* cm"'^ denotes the transition to self-shielding. Ha- 
los with M > 10® Mq begin to self-shield Lyman- Werner radiation, 
where most of the LW radiation is self-shielded in Tyir = 10* K 
halos. 



3.1. Number of Star Forming Halos 

Gas in halos with masses <10^ Mq is evacuated by 
~30 km s^^ D-type front, leavin g a diffuse (1 cm~'^) 



and war m (3 x 10"^ Kl medium (iWhalen e t al.' '2004' 
Kitayama et al.l |2004 lYoshida et all 120071 : lAbcl ct al. 
2007t) . The aftermath of SN explosions in relic H II 
regions was explored in spherical symmetry in one - 
dimensional calculations bv lKitavama fc Yoshidal (|2005D . 
Even without SNe, star formation is suppressed for ^^100 
Myr before gas is reincorporated into the potential well 
of this early galactic progenitor. PISNe provide an extra 
^ 10^^ erg of thermal energy and can evacuate halos up 
to 10^ Mq. Hence star formation within these low-mass 
halos are highly dependent on their sta r format ion and 
merger histories, as illustrated by Yoshi da"et al.l (j20G7). 

Table [2] summarizes the global properties of the most 
massive halo the time of collapse in the HHe calculations. 
Approximately 5-20 Pop HI stars form in the progeni- 
tors, whose original gas structures are nearly destroyed 
by radiative feedback, of the lO"* K halo. More specifi- 
cally at z = 15.9 (16.8), there are 19, 29, and 24 stars 
that form after redshift 30 in the SimA-RT, SimB-RT, 
and SimB-SN runs, respectively. The most massive halo 



Name 


N*(<rvir) 


N*(<3rvir) 


Mgas/Mtot 




SimA-HHe 






0.14 


0.010 


SimA-RT 


11 


12 


0.064 


0.030 


SimB-HHe 






0.14 


0.010 


SimB-RT 


22 


26 


0.089 


0.014 


SimB-SN 


4 


10 


0.046 


0.038 



Note. — Col. (1): Simulation name. Col. (2): Number 
of stars hosted in the halo and its progenitors. Col. (3): 
Number of stars formed in the Lagrangian volume within 
3rvir of the most massive halo. Col. (4): Baryon fraction 
within rvir. Col. (5): Baryonic spin parameter [Eq. [5]. 



and its progenitors have hosted 11, 22, and 4 stars in 
the same simulations. When we look at the Lagrangian 
volume contained within three times the virial radius of 
the most massive halo, there have been 12, 26, and 10 
instan ces of star formation. For comparison, iGreif et al.l 
(|200l find an an upper limit of 10 Population HI stars 
that are hosted in the progenitors of an early galaxy sim- 
ilar to the ones presented in this work. 

3.2. Global Nature of Objects 

In addition to providing the first ionizing photons and 
metals to the universe. Pop HI stars change the global 
gas dynamics of Tvir '~ 10'' K star forming halos. Figures 
|2]and|3|compare the structure of the most massive halo in 
all simulations, depicting density-squared weighted pro- 
jections of gas density and temperature. All of the halos 
have a virial mass of 3.5 x 1O^M0. The models with 
neither star formation nor H2 chemistry show a centrally 
concentrated, condensing Tvir = 10^ K halo with its asso- 
ciated virial heating. In comparison, feedback from pri- 
mordial star formation expels the majority of the gas in 
low-mass star forming progenitors. The combination of 
these outflows and accretion from filaments induces the 
formation of an inhomogeneous medium, where the ra- 
diation anisotropically propagates, creating champagne 
flows in the directions with lower column densities. The 
temperature projections illustrate both the ultraviolet 
heated 10"* K) and optically thick, cool 10^ K) 
regions in the host halo and IGM. Also note the nearby 
substructure is photo-evaporated by the stars hosted in 
the most massive progenitor. SN explosions alter the gas 
structure by further stirring and ejecting material after 
main sequence. Next we quantify these visual features 
with the baryon fraction of the halo and an inspection of 
phase diagrams of density and temperature. 

Because the gas in the progenitors is mostly evacuated, 
the baryon fraction of high-redshift star forming halos are 
greatly reduced. In halos with masses < 10^ M©, the 
baryon fraction lowers to 5 x 10~'^ ten million years afte r 
the star's death without a SN (cf. lYoshida et all 120071) . 
When we include SNe, the baryon fraction decreases fur- 
the r to 1 X 10~^ six million years after the explosion 
(cf. iKitavama fc Yoshidall2005l) . This is in stark contrast 
with the cosmic fraction Q.i,/U,m = 0.17. 

Within these shallow potential wells, outflows from 
stellar feedback impede subsequent star formation until 
sufficient gas is reincorporated, occurring through merg- 
ers and smooth IGM accretion. After the halo mass 
surpasses ~3 x 1O^M0, total evacuation does not oc- 
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z = 15.9 
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Fig. 2. — Density-squared weighted projections of gas density 
{left) and temperature (right) of the most massive halo in simu- 
lation A. The field of view is 1.2 proper kpc. The top row shows 
the model without star formation and only atomic hydrogen and 
helium cooling. The bottom row shows the same halo affected by 
primordial star formation. Note the filamentary density structures, 
clumpy interstellar medium, and the counter-intuitive effect that 
feedback leads to lower temperatures. 



cur but significant outflows are still generated. Near 
the same mass scale, multiple sites of star formation oc- 
cur in the same halo. In our simulations, these stars 
rarely shine simultaneously since their lifetimes is sig- 
nificantly shorter than the Hubble timescale on which 
the halos are assembled. We neglect H2 self-shielding, 
which may affect the timing of star for mation but un- 
likely not the global star formation rate (|Ahn fc Shapirol 
|2007() . Here the nearby overdensity survives the blast 
of UV photo-heating from the previous star and con- 
denses to form a star a few Myr afterwards. This sce- 
nario of adj acent star f o rmati on is similar to the one 
presented in lAbel et all (|2007t l. but the multiple sites 
of star formation are caused by H2 and Lyg cooling 
in central protogalactic shocks (Shapiro fc Kandll987f ) 
and aided by addit i onal H2 cooling in ionization fronts 
(iRicotti et all 1200 ll : lAhn fc Shapirol l2007t IWhalen et all 
|2008|) . not from residual cores from a recent major 
merger. 

After the first star forms in the most massive progen- 
itor, the halo slowly regains gas mass mainly through 
merger events of halos that have not experienced star 
formation. We show the evolution of the baryon fraction 
of simulation B in the top section of Figure H) In SimB- 
RT, ionization fronts created from stars in halos with 
Mvir 3 X IO^AIq cannot expel the majority of baryons 
in the halo but only generate outfiows. In SimB-SN, sim- 
ilar recovery occurs but the additional energy from the 
three SNe at redshift 20 evacuates the halo once again to 
a baryon fraction of 0.02. When the most massive halo 
reaches Tvir ^ 10** K, the baryon fraction within the virial 
radius has only partially recovered to 0.064, 0.089, and 
0.046 in the SimA-RT, SiniB-RT, and SimB-SN calcula- 




FlG. 3. — Same as Figure[2]for simulation B. Here the bottom row 
shows the halo with primordial stellar feedback and supernovae. 



tions. Without any stellar feedback, these fractions are 
0.14 in the HHe runs. 

During central gaseous collapses in halos, the gravita- 
tional potential in the inner '^50 pc becomes baryon dom- 




25 

Redsiiift 

Fig. 4. — The top graph shows the total mass (thick black line) 
of the most massive progenitor in simulation B. The thin lines 
depict the baryon fraction, A/gas /Aftot, of simulation B without 
star formation (blue solid), with star formation (red dashed), and 
with SNe (green dot-dashed). The bottom graph shows the spin 
parameter (eq. |5ll of the same halo in the dark matter (thick line) 
and gas (thin lines, same legend as top graph). 
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Fig. 5. — Radial profiles of DM density (top panel) and enclosed 
DM mass (bottom panel) of the most massive halo in simulation B 
at z = 16.8 with an adiabatic equation of state (solid), no star for- 
mation (dotted), star formation and feedback (dashed), and plus 
SN feedback (dash-dotted). The reduced central baryon overden- 
sities in the progenitors lead to lower central DM densities. 
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Fig. 6. — Mass-weighted radial profiles of density (top) and 
temperature (bottom), centered on the densest DM particle. The 
columns show data from SimA-RT (left), SimB-RT (middle), and 
SimB-SN (right). Note how the cool and warm gas phases coexist 
at similar radii throughout the halo. 



inated where the baryon density is greater than the DM 
density, leading to a co n traction of the DM inner halo 
(|Blumenthal et al.lll986t iGnedin et a"l]|2004D . If stellar 
feedback evaporates this overdensity, the central poten- 
tial will not be as deep during the assembly of early dwarf 
galaxies, possibly resulting in the DM being not as cen- 
trally concentrated. We plot radial profiles of DM den- 
sity and enclosed DM mass of the most massive halo in 
simulation B at z = 16.8 in Figure E] The dark matter 
density is decreased in the inner 20 pc (O.OSrvir) up to 
a factor of 5 with stellar feedback alone. The effect is 
exacerbated by the additional feedback from SN explo- 
sions inside 50 pc (O.lrvir), and the central DM densities 
decrease by another factor of 2, similar to adiabatic cal- 
culations. 

These outflows also create inhomogeneities in and 
around halos and increase the baryonic spin parameter 



where Lg, Eg, and Afg are the total baryonic ang ular mo- 
mentum, kinetic energy, and mass of the system ("Peeblesl 
[l971). This is basically the ratio of rotational to gravi- 
tational energy of the system. The DM spin parameter 
A uses the total DM angular momentum, kinetic energy, 
and mass of the system. 

Interestingly the halo experiences significant fluctua- 
tions in Xg when significant outflows are generated by 
stellar feedback. We plot its evolution in the most mas- 
sive progenitor in the bottom section of Figure S) With- 
out feedback (HHe) after redshift 25, the gas and DM 
spin parameters are approximately equal and follow the 
same trends. However when star formation (SimB-RT) 
is included, the spin parameter is increased by an order 
of magnitude after the first star. It then decays over the 
next 40 Myr but increases again after the second star 
in the halo at redshift 26. Xg continues to be up to a 
factor of two higher than without star formation after 
the star can no longer expel most of the gas from the 
halo. In SimB-SN, these effects are even more apparent, 
especially at redshift 20 during an episode of three SNe 



evacuating the most of the gas from the halo. 

At the time of collapse in the HHe runs, A = 
0.022 (0.013) and is marginally lower than the av- 
erage (A) 0.04 found in numerical simulat ions 

(|Barne s fc Efstathiou 1987^. 'Eise nstein fc Locb 19951). It 
is not affected by stellar feedback as DM dominates the 
potential well. Without star formation, the baryonic spin 
parameter Ag = 0.010 (0.010) and is slightly lower than 
A. However with stellar and SNe feedback, Ag increases 
up to a factor of 4. The effect is smaller without SNe but 
still significant, raising Ag to 0.030 (0.014). 

The increase of Ag could be caused by the forces gen- 
erate d in the H II region and SNe blastwaves ()Abel et al.l 
l2001h . Because these events can expel gas from the po- 
tential well and (A) ~ 0.04, these forces can be up to 
25 times greater than the cosmological tidal torques usu- 
ally associated with the angular momentum of galaxies 
(|Peebleslll969D . The angular momentum associated with 
ionization fronts and blastwaves would have to be al- 
most opposite of angular momentum vector of the halo 
to slow the rotation. Therefore we expect these outflows 
to produce on average an overall increase the angular 
momentum of galaxies. In the assembly of early dwarf 
galaxies, a fraction of the expelled gas experiencing large 
scale torques longer falls back, now with a higher specific 
angular momentum, to the galaxy and increasing its spin 
parameter further. 

3.3. ISM Phases 

The combination of molecular cooling, stellar feedback, 
and SN explosions create a multi-phase ISM in star form- 
ing halos. These phases are interspersed throughout the 
halo. They are marginally seen in the temperature pro- 
jections in Figures [5] and [31 However they are better 
demonstrated by the mass- weighted radial profiles in Fig- 
ure [B] and density-temperature phase diagrams in Fig- 
ure [T] The radial profiles are centered on the densest 
DM particle. For a given radius within the halo, the gas 
density can span up to 6 orders of magnitude, and the 
temperature ranges from '--'100 K in the cool phase to 
30,000 K in the warm, ionized phase. Below we describe 
the different ISM phases at redshift 15.9 and 16.8 for 
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Fig. 7. — Mass-weighted p-T phase diagrams of a sphere with 
radius 1 kpc, centered on the most massive halo in SimA-RT (top), 
SimB-RT (middle), and SimB-SN (bottom). At T > 10* K, one 
can see the H II regions created by current star formation. The 
warm, low density (p < 10~^ cm~'^) gas in SimB-SN are contained 
in SNe shells. 



simulation A and B, respectively. 

Cool phase — The relatively dense (p > 100 cm^'^) 
and cool (T > 1000 K) gas has started to condense 
by H2 cooling. Current star formation dissociates H2 in 
nearby condensations through LW radiation in our sim- 
ulations. In many cases, especially when Mvir ^ 10'' Mq, 
nearby clumps remain cool and optically thick. After the 
star dies, H2 formation can proceed again to form a star 
in these clumps. There are two other sources of cool gas. 
First, the filaments are largely shielded from being photo- 
heated and provide the galaxy with cold accretion flows. 
Second, after the main sequence, the material within the 
expanding shell either from a D-type front or SN blast- 
wave cools through adiabatic expansion and Compton 
cooling to temperatures as low as 100 K, which is seen 
in the p-T phase diagram at very low densities. 

Warm, neutral phase — Gas that cools by atomic hy- 
drogen line transitions, but not molecular, has T ^ 8000 
K and densities ranging from lO^'^ to 10^ cm^^. Gas in 
relic H II regions and virially shock-heated gas compose 
this phase. 

Warm, ionized phase — In the SimB-RT simulation at 
the final redshift, a Pop III star is creating an H II region 
with temperatures up to 30,000 K. It is 20 pc from the 
halo center. The H II region does not breakout from this 




log p [cm ] 



Fig. 8. — The same phase diagram of SimB-SN in Figure [7] 
but colored by mean metallicity and scaled to show the hot, X-ray 
phase. These data are 45 kyr after a pair-instability SN. The SN 
remnants that are warm and diffuse (p < 10"'^ cm~^) have solar 
metallicities or greater. The majority of the ISM has metallicities 
^j^Q— 2.5 sQiar. The densest, collapsing material has a metallicity 

T^ir ^ 10* K halo. 

Hot, X-ray phase — The 3 x 10^^ ergs of energy de- 
posited by one PISN in the SimB-SN simulation heats 
the gas to over 10* K immediately after the explosion. 
Figure [8] shows the density and temperature of the ISM 
45 kyr after a SN, where the adiabat of the hot phase is 
clearly visible at T > 10^ K. A blastwave initially travel- 
ling at 4000 km s^^ sweeps through the ambient medium 
during the free expansion phase. The gas behind the 
shock cools adiabatically and through Compton cooling 
as the SN shell expands. 

IWise fc Abell ()2007[ ) investigated the generation of tur- 
bulence during virialization without stellar feedback. 
The supersonic nature of turbulence in these halos re- 
mains even with Pop III stellar feedback and H2 cool- 
ing. In Figure [9l we plot the turbulent Mach number, 
'Vrrns/cs, in the most massive halo. Here Vrms is the 
three-dimensional rms velocity relative to the mean ve- 
locity of each spherical shell, and is the sound speed. 
Turbulent Mach numbers reach 6 in SimA-RT but is only 
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Fig. 9. — Turbulent Mach numbers for SimA-RT (solid), SimB- 
RT (dashed), and SimB-SN (dotted). The Mach numbers in sim- 
ulation B are lower than simulation A because of the higher tem- 
peratures created by stellar feedback that occurred shortly before 
the displayed data. 
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Fig. 10. — Density-squared weighted projections of metallicity 
{color) and gas density {black and white) for SimB-SNe at 2; = 16.8 
of the inner 8.6 proper kpc of the simulation. 



1-2 in SimB-RT and SimB-SN, where the recent stel- 
lar feedback in the halo has photo-heated the gas. In 
contrast, sufficient time has elapsed since the previous 
episode of star formation in SiniA-RT, thus allowing the 
gas to cool by H2 to ~300 K. Clearly these values are 
dependent on the stellar feedback within the halo and 
possibly its merger history, hence this range may be rep- 
resentative of turbulent Mach numbers one would find in 
early dwarf galaxies. 

3.4. Metallicity 

Metallicities of second and later generations of stars 
depend on the location of previous SNe. Figure [TOl shows 
a projection of metallicity for the inner 8.7 proper kpc of 
the SimB-SNe at z = 16.8, which is overplotted against 
the density-squared weighted projection of gas density. 
We also show the metallicity projection for the region 
surrounding the most massive halo in Figure 1111 This 
projection depicts the data in a cube with a side of 1.2 
proper kpc, centered on the halo. The contours in this 
figure mark the number density of gas. In the SimB-SN 
calculation, outflows carry most of the SN ejecta to radii 
up to '--^l proper kpc after 30 Myr. Interestingly they 
approximately fill the relic H II region and expand little 
beyond that. The low density IGM marginally resists the 
outflows, and it is preferentially enriched instead of the 
surrounding filaments and halos. 

It should be noted that this calculation is an upper 
limit of metallicities since all stars end with a PISN. The 
mixing and transport of the first metals is a fundamental 
element of the transition to Pop II stars and is beneficial 
to study in detail. All metallicities quoted are in units 
of solar metallicity. The metallicities also scale approxi- 
mately linearly with metal yield of each SN because we 
treat the metal field as a tracer field that is advected 
with the fluid flow. We quote the metallicities according 
to this scaling. 




Fig. 11. — Density-squared weighted projections of metallicity, 
centered on the most massive halo in SimB-SNe at z = 16.8. The 
contours depict the number density of baryons for n = (4, 13, 50, 



320) 



The field of view is 1.2 proper kpc and the same as 



Figure|3] The projection shows data in a slab that are 1.2 prope 
kpc thick. 



We plot the mean metallicity of the most massive pro- 
gentior of this dwarf galaxy in Figure [T^] as a function 
of redshift. The first star enriches the gas depleted host 
halo, and the metallicity decreases as metal-free gas is 
incorporated through mergers and IGM accretion. The 
temporary increase in metallicity at z = 20 is associated 
with the three stars that form in succession. 

When the most massive halo reaches a virial tempera- 
ture of 10** K at z = 16.8, metals are thoroughly mixed 
in the halo, and its mean metallicity is lO^^'^/go, where 
/so is in units of 80 Mq of metal ejecta per PISN. Turbu- 
lence created from the dynamic assembly of this Tvir = 




26 24 22 

Redshift 

Fig. 12. — Mean metallicity (solid) and baryon fraction (dashed) 
in the most massive progenitor in SimB-SN. The first star enriches 
the diffuse gas above solar metallicity, which then decreases as pris- 
tine gas falls into the halo through accretion and mergers. Notice 
that the metallicity stays roughly constant at 10~^ when no star 
formation is occurring. 
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Fig. 13. — Mass-weighted p-T phase diagrams of a sphere with radius 1 kpc, centered on the most massive halo in SimB-SN at redshift 
21 (left) and 16.8 (right). The halo mass at z = (21, 16.8) is 1.1 X lO'' Mq and 2.9 X IO'^Mq, respectively. The pixel intensity is determined 
by the average of the ratio f of the cooling rate with metal-line cooling and the cooling rate with primordial chemistry. At z = 21, we 
underestimate cooling rates in gas with T> 1000 K and p ~ 0.1 cm"'' up to a factor of 10, which only composes ~2% of the halo gas 
mass. However, as the gas becomes more enriched, metal-line cooling becomes dominant over H, He, and H2 cooling in cool gas and old 
SN remnants (T > 10'' K, p < 10~^ cm~^). However, the latter has a cooling time greater than the Hubble time. 



10^ K halo that involves outflows generated by feedback, 
cold accretion through filaments, halo mergers, and viri- 
aliza tion (Paper I) app ears to thoroughly mix the met- 
als (jPan fc Scald l2007l ). As stars form in this halo, the 
metallicity of this halo fluctuates around this value be- 
cause stars continue forming but ejecting most of their 
metals into the IGM. Also the filaments are still mostly 
pristine and provide a source of nearly metal-free cold 
gas. The volume averaged metallicity of the enriched 
IGM ((5 < 10) is 10~^ °/8o, compared to the filaments 
and halos (6 > 10) that ar e less enriched with 10~^-^/8o 
(cf. iCen fc Riauelmell2008D . Considering the total mass 
of heavy elements, there are 360 Mq of heavy elements in 
the IGM, about 40% of the total ejecta from all PISNe. 
The metal volume filling fraction is 3.5% of the volume 
where we allow star formation to occur. This percentage 
should be higher than the cosmic mean at this redshift 
because this comoving volume of (205 kpc)'^ is biased 
with an overdensity S = p/ p = 1.8. Thus star formation 
rates are greater than the mean since there are more 
high-cr peaks, and the metal filling fraction should scale 
with this bias. 

The p-T phase diagram in Figure [5] is colored by the 
mean metallicity of the gas. These data are taken imme- 



diately before the formation of a star with a metallicity 
of 10~"^'^/8o. There are three distinct metallicity states 
in the halo. The majority of the gas in the halo has a 
density between 10"'^ and 1 cm~'^. This gas has a mean 
metallicity of 10~^'^/8o- At higher densities, the metal- 
licity is slightly lower at 10~^'^/8o. The same preferential 
enrichment of diffuse regions may have caused the lower 
metallicities in this dense cloud. The third phase is the 
warm, low density {p < 10^^ cm~'^) gas that exists in re- 
cent SN remnants and has solar metallicities and greater. 

The hot phase produced by a SN is super-solar. The 
high-density tail of the ejecta is the SN shell and is 
mixing with the lower metallicity ambient medium. As 
the ejecta expands and cools, it will contribute to the 
warm, low density, solar metallicity material, whose cool- 
ing time is greater than the Hubble time. 

Although we do not consider metal-line cooling in the 
enriched material, it is useful to determine the impor- 
tance of metal-line cooling in early galaxies. As stated 
above, heavy elements within the most massive halo are 
well-mixed by turbulence and should contribute to the 
cooling of the ga s as it is enriched above ^1 0"^^^ (e.g. 
iMaio et al.ll2007t ISmith fc Sigurdssonl 120081 ). We calcu- 
late cooling rates with contributions from metal-line cool- 
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Fig. 14. — Left: Fraction of mass (black solid) and volume (red dashed) of the sphere with a radius of 1 kpc, centered on the most 
massive halo at 2 = 21 (left) and 16.8 (right), that would have a high cooling rate by a factor of if we would consider metal-line cooling. 
Right: Mass-weighted (top) and volume-weighted (bottom) histograms of cooling rates of the same volume at redshift 21 (left) and 16.8 
(right). The solid black lines are the cooling rates with primordial abundances, whereas the red dashed lines consider metal-line cooling. 
The increased cooling at 10"'^^ erg s~^ is caused by metal-enriched cool gas, and the increased cooling above 10~^^ originates from warm 
( Si 10^ K) gas, mainly in SN remnants. 



ing in the most massive halo in SimB-SN. We stress that 
the following cooling rates from metals are upper lim- 
its since every Population III star produces a PISN. For 
gas with temperatures below 10* K, we calculate the ad- 
ditional cooling from carbon, oxygen, and silicon using 
the rates and reactions from Glover & Jappsen (2007). 
For hotter gas, we u se th e metal cooling rates from 
[Sutherland fc Dopital (| 19931) . i.e. AA(Z) = Anet(2') - 
Anetf-^ = 0), where A net is the net cooling rate listed in 
[Sutherland & Dopita' and Z is the metallicity relative to 
solar. Figure [T3l compares the cooling rates with only 
primordial abundances (i.e. lAbel et anil997f ) and with 
metal-line cooling in a density-temperature phase plot. 
We plot the data from a sphere of 1 proper kpc, cen- 
tered on the most massive halo when it forms its second 
star {z — 21) and the final redshift (z = 16.8). The 
pixel intensity represents the ratio C — Amotais/Ao, the 
cooling rate with metal-line cooling and with primordial 
abundances. The gas that experiences the greatest metal 
cooling compared to cooling from H, He, and H2 is either 
in old SN remnants and H II regions that have T > 10* 
K or in cool (T < 1000 K) and diffuse (p - 1 cm'^) gas 
where fine-structure metal lines start to dominate over 
II2 cooling. The warm (T « 10* K) ISM cools by atomic 
hydrogen line cooling and metal line cooling and the re- 
cent SN remnant {T > 10^ K) cools mainly by Compton 
cooling. 

We plot the cumulative sum of mass and volume that 
as a function of C at z = (21, 16.8) in the left panel 
of Figure [TH Approximately 12% (41%) of the mass 
in the most massive halo at redshift 16.8 (21) have no 
significant metal-line cooling. At the final redshift, over 
a quarter of the gas mass has cooling rates increased 
over an order of magnitude. However z = 21, this mass 
fraction drops to 4%. In the right panel, we also plot 
the probability density function of cooling rates in both 
the enriched and primordial cases. We see the features 
from H II regions, SN remnants, and cool gas in Figure 



[T5I also in Figure [TH The shift of the first bump that 
was originally centered at 10~^^ erg was caused by 
the cool metal-enriched gas, and the shift at higher rates 
occurs in H II regions. Taken as an upper limit of metal 
enrichment, it appears that metal cooling considerably 
alters the cooling characteristics of the gas within these 
early dwarf galaxies and should be accounted for in future 
simulations that study the early dwarf galaxy formation 
or the transition to Pop II star formation. 

As the metallicities in early dwarf galaxies increase, 
dust absorption may become important. The dust ex- 
tinction cross-section in the Milky Way is ad,MW — 3 x 
10~^ ^ cm^ per hydrogen nucleus (jHollenbach fc McKed 
[T97I ICardelh et al.l[l98l . If we assume the dust ex- 
tinction properties do not change at lower metallicities, 
the cross-section scales with metallicity, giving — 
{Z/ZQ)ad,MW ■ The optical depth to dust absportion 
is simply 

f J, fd.MW ^metals /„n 

Td= riHOd al = — ■ , (6) 

J ^0 

where A^motais is the column density of metals and Zq — 
0.0204. In Figure [T51 we calculate the distribution of 
A'^mctais between the center of the most massive halo and 
a sphere of radius 3 kpc with our ray tracer, exactly like 
in the H2 case in Figure [TJ Column densities of metals 
are distributed between 2 — 5 x lO^'* cm^, corresponding 
to an upper limit of Td = 0.3 — 0.7. Hence dust shielding 
should not be significant in the formation of early dwarf 
galaxies. More massive halos, however, may be subject 
to dust shielding as the metallicity and total hydrogen 
column density of the halo increase. 

3.5. Metal-enriched Star Formation 

The star formation times of SimB-SN is depicted in 
Figure [TH by plotting the total mass of the host halo ver- 
sus the formation redshift. The different symbols repre- 
sent the metallicity of the star. Before redshift 20, zero 
metallicity stars form in low-mass halos, whose masses 
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Fig. 15. — Distribution of metal column densities between the 
center of the most massive halo and a sphere of radius 3 proper 
kpc at 2 = 16.8 in SimB-SN. 
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Fig. 16. — Star formation times of SimB-SN. The y-axis shows 
the total mass of the host halo. The symbols indicate the different 
mctallicities of the stars and are labeled in the legend. 



increase from 5 X W^Mq to 2 x IO^Mq due negative 
feedback from photoevaporation of low-mass halos. The 
amount of photoevaporation in nearby halos is appar- 
ent in the density projections in Figures [2H3] where most 
nearby substructure is lacking in the simulations with 
star formation. This will further suppress star forma- 
tion already hindered by the LW radiation background 
(jMachacek et al.ll2001[ ). 

The first instance of a metal enriched ([Z/H] > -6) 
star occurs at z = 20.7 with a metallicity of 10~*'°/8o in 
the most massive halo that has a mass of 9.5 x IO^Mq. 
This star was trigge red by a SN b lastwave, not in the 
IGM as envisaged bv lFerraral ()1998[ ) but within the same 
halo. This SN explodes 470 kyr before the star and pro- 
vided the majority of heavy elements for the formation 
of this enriched star. The SN of the second star in the 
halo triggers yet another round of star formation with 
a metallicity of lO^'^'^/so only 80 kyr afterwards. The 
aggregate energy from these three SNe expel most of the 
gas from the potential well. The most massive halo does 
not form any stars until z = 16.9 at which it undergoes 
two episodes of star formation both with a metallicity of 

10-2-780. 

These metal-enriched stars could exceed the 

critic al metallicity of ~ 10 "'* Zr^ (iBromm et ahl 

20011: ISchneid cr ct al. 2006; iJappsen et all l2007bl : 



Smith fc Sigurdsson 2007) . If this were to happen, the 
metal-enriched g as would collapse and fragment into 
a stellar cluster (jClark et al.ll2008h . However it is too 
preliminary to determine the differences in radiative and 
chemical feedback between a Population HI star and a 
Population II stellar cluster. As a point of reference, a 
lO^"^ Zq star cluster with a Salpeter IMF would need 
to host 3500 Mq of stars to produce equal amounts of 
hydrogen ionizing photons of what alTOMQ Population 
HI star produces (jSchaereil l2002l |2003| ). Assuming a 
typical SN explosion energy of 7 x lO''* erg Mg^ in a 
stellar cluster with a Salpeter IMF, one would similarily 
need 3900 Mq of stars to equal the explosion energy of 
a 170 Mq metal-free star. 

It should be noted that four stars in our calculation 
form with a trace of metals ([Z/H] ^ -6). Three in- 
stances happen when a SN blastwave overtakes a nearby 



halo. Here the blastwave leaves the high density mate- 
rial fairly pristine and the surrounding IGM is enriched 
(|Cen fc Riouelmd 12001 ■ 

In three halos with masses 1 — 4 x IO^Mq, star forma- 
tion is triggered shortly (< 3 Myr) after the death of a 
previous star in the same halo. The second star in the 
halo forms in a density enhancement that is caused by an 
ioniza tion front instability (e.g., see IWhalen fc NormanI 
l2008airbh when a SN blastwave overtakes it. The size of 
the overdensity is small enough (~30 pc) so the heavy 
elements can mix into the high density material within 
a free-fall time. In these three halos, this results in stars 
with metallicities of [-5.3, -6.4, -4.3] (Myieid/Mg) that 
form at z = [19.9, 19.1, 17.8], respectively. These metal- 
licities are uncertain due to the timing of star forma- 
tion since we neglect H2 self-shielding. If the density 
enhancement was shielded from the LW radiation from 
the neighboring star and H2 cooling was still efficient, 
the star could have formed before the SN blastwave en- 
riched the cloud. In that case, the resulting star would 
have been metal-free. Hence the metallicities of these 
induced star forming regions should be considered with 
caution. Studies with a more accurate treatment of the 
H2 self-shielding could better capture the temporal se- 
quence of metal enrichment and ionization front insta- 
bilities. Given the uncertainties in collapse times from 
our simple star formation algorithm, masses, and feed- 
back parameters, it is improbable that a better H2 line 
transfer calculation would lead to significantly realistic 
results. 

4. DISCUSSION 

We find that the combination of Pop HI stellar feed- 
back and continued H2 cooling in T^n- < lO'' K halos 
alters the landscape of high-redshift galaxy formation. 
The most drastic changes are as follows: 

1. Dynamic assembly of dwarf galaxies — A striking 
difference when we include Pop HI radiative feedback 
are the outflows and gas inhomogeneities in the halos 
and surrounding IGM. The outflows enrich the IGM and 
reduce the baryon fraction of the 10** K halo as low as 
0.0 5, much lower than the cosmic fraction ^Ib/flM — 0.17 
(cf. lYoshida et al.ll2007l : [Abel et al.ll2007t ). This substan- 
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tially differs from the current theories of galaxy forma- 
tion where relaxed isothermal gas halos hierarchically 
assemble a dwarf galaxy. For instance in simulation A, 
there are remarkable filamentary structures and a clumpy 
ISM. Furthermore, Pop III feedback increases the total 
baryonic angular momentum of the system up to a factor 
of 3 without SNe and up to 5 with SNe. 

2. Pop III sphere of influence — Pop 111 feedback is 
mainly a local phenomenon except its contribution to 
the UVB. How far its H II region, outflows, and metal 
ejecta (if any) will predominately determine the charac- 
teristics of the next generation of stars. Highly biased 
(clustered) regions are significantly affected by Pop HI 
feedback. The first galaxies will form in these biased re- 
gions and thus should be significantly influenced by its 
progenitors. 

3. Dependence on star forming progenitors — Al- 
though our calculations with SNe only provided an upper 
limit of metal enrichment, it is clear that the metallicity, 
therefore metal-line and dust cooling and metal-enriched 
(Pop 11) star formation, depends on the nature of the 
progenitors of the dwarf galaxy. If the galaxy was as- 
sembled by smaller halos that hosted a Pop HI star that 
did not produce a SN, e.g., the galaxy would continue to 
have a top-heavy initial mass function (IMF). 

4. Complex protogalactic ISM — The interplay be- 
tween stellar and SNe feedback, cold inflows, and molecu- 
lar cooling produce a truly multi-phase ISM that is remi- 
niscent of local galaxies. The cool, warm, and hot phases 
are interspersed throughout the dwarf galaxy, whose tem- 
peratures and densities can span up to 6 orders of mag- 
nitude at a given radius. 

5. Metallicity floor — When the halo is massive enough 
to host multiple sites of star formation, the metal ejecta 
does not significantly increase the mean metallicity of 
the host halo. There seems to be a balance between 
(a) galactic outflows produced from SNe, (b) inflowing 
metal-enriched gas, (c) inflowing pristine gas, and (d) 
SNe ejecta that is not blown out of the system. In our 
high yield models, the metallicity interestingly fluctuates 
around Zq in the most massive halo when this bal- 
ance occurs at and above mass scales ~10^Mq. 

Clearly the first and smallest galaxies are complex en- 
tities, contrary to their low mass and generally assumed 
simplicity. Our calculations refiect the important role of 
Pop III stellar feedback in early galaxy formation. 

These high-redshift galaxies have a ~5-15% chance of 
being undisturbed by merger s until the present day, be - 
ing "fossils" of reionization (|Gnedin fc Kravtsovl |2006| ) . 
Dwarf spheroidals (dSph) galaxies are some of the dark- 
est galaxies in the universe, ha ving high mass-to- light 
ratios up to 100 (jMateol Il998f) . Gas loss in dSph's 
close to the Milky Way or M31 can be expl ained by gas 
tidal stripping during orbital encounters (jMaver et al.l 
l2007f) . However there are some galaxies (e.g. Tucana, 
Cetus) removed from both the Milky Way and An- 
dromeda galaxies and cannot be explained by tidal strip- 
ping. In addition to ultraviolet heating from reionization 
(|Bullock et al.ll2000l-[Su"sa fc Umemurall2004D a nd intrin- 
sic star formation (|Mac Low fc Ferraral Il999l ). perhaps 
stellar feedback from Pop III stars influenced the gas- 
poor nature of dSph's. Even at the onset of widespread 
star formation in the objects studied here, the baryon 
fraction can be three times lower than the cosmic mean. 



and the dwarf galaxy may never fully recover from the 
early mass loss. This initial deficit may play an impor- 
tant role in future star formation within these low-mass 
galaxies and could help explain the large mass to light 
ratio in isolated dSph's. 

With the radiative and chemical feedback from the pro- 
genitors of the early dwarf galaxies, we have an adequate 
set of cosmological "initial conditions" to study the tran- 
sition from Pop III to Pop II stars. In this setup, the 
current metal tracer field can be used to include metal 
line cooling. Dust cooling may induce fragmentation of 
solar mass frag ments at metallicities as low as ~10~^ at 
high densities (jSchneider et al.|[2006f) . However, metal- 
line cooling might n ot be important at the se low metal- 
licities in diffuse gas. iJappsen et al.l (|2007a D showed that 
metal-line cooling at metallicities below 10~^Zq in low 
density gas does not significantly affect t he dynamics of 
a collap sing halo. In a companion paper, iJappsen et al.l 
(|2007bf) found that the fragmentation of a metal-poor 
{Z = 1Q~^-^Zq) collapsing object may depend more on 
the conditions, e.g. turbulence and angular momentum 
distributions, created during the assembly of such a halo 
than some critical metallicity. 

Perhaps when the protogalactic gas cloud starts to host 
multiple sites of star formation, the associated SNe pro- 
duce sufficient dust in order for a transition to Pop II. 
In lower mass halos, the SN ejecta is blown out of the 
halo, and future star formation cannot occur until ad- 
ditional gas is reincorporated into the halo. However in 
these halos with masses 10 ''M©, the SN does not to- 
tally disrupt the halo. A fraction of the SN ejecta is con- 
tained within halo and could contribute to subsequent 
sites of star formation. SN ejecta and associated dust 
could instigate the birth of the first Pop II stars. 

The metallicities of the first generation of metal- 
enriched stars could depend on the metal mixing 
timescales in these early dwarf galaxies. Fortunately 
the dispersion of heavy elements from SNe into the 
Galactic ISM is a rich fi eld of study (for a review, see 
iScalo fc Elmegreenll200l . MetaUicity dispersions in lo- 
cal stellar clusters show fluctuations of 5-20% around 
the mean, which suggest t hat the ISM is well-mixed (e.g . 
Edvardsson et al.l 119931 : iGarnett fc Kobulnickvl 120001 : 
Reddv et al.]l2003l ). In addition, very metal-poor ([Fe/H] 
< -2.7) halo stars have very little scatter in elemental 
abundances, suggesting that the dispersal of the first 
metals that formed low mass stars originate d from single 
starb ursts instead of single SN explosions (jCavrel et al.l 
I2OOI . 

In the ISM, laminar flows and turbulence driven 
by SN explosions provides the impetus of metal mix- 
ing on large scales, which then cascades to smaller 
and smaller length scales eventually reaching a length 
scale associated with a Reyn olds number Re — 1 (e.g. 
Ide Avillez fc Mac Lowl 120021 ). In structure formation, 
turbulence can also arise durin g the virialization of cos- 
mological halos (iWise fc Abefl 120071 : iGreif et al.l l20?)8^. 
After the turbulent cascade creates metallicity gradients 
on small enough scales, mol ecular diffusion will ho mog- 
enize the heavy elements. Ide Avillez fc Mac Low! per- 
formed AMR simulations that include turbulent diffu- 
sivity to study the mixing timescales in the ISM. In their 
simulations and the ones presented here, numerical dif- 
fusion provides the majority of the mixing at the res- 
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olution limit. They find that the mixing time is fairly 
independent of the diffusion scale and could depend on 
some inertial scale. The origin and magnitude of diffu- 
sion does not significantly affect the gas on large scales. 
In their resolution study, mixing timescales only decrease 
by 20% when the resolution is increased by a factor of 4. 
Thus we could be overestimating the mixing timescales 
up to a factor of a few because it is difficult to resolve 
the smallest turbulent scale in cosmological simulations. 
Furthermore, we expect metals that are newly incorpo- 
rated into the galaxy to be stirred by the turbulence in 
assembling halos, which is sustained in virialization and 
major mergers on the largest scales. 

As discussed above, the metallicity of the most mas- 
sive halo fluctuates around IO^^Zq. This is intriguingly 
the same value as a sharp cutoff in stellar metallicities 
in fo ur local dSph's: Sculptor, Sextans, For nax, and Ca- 
rina (|Tolstov et al.ll2004HHelmi et al.ll2006D . This differs 
with the galactic halo st ars, whose metal-poor ta il ex- 
tends to ZIZq = 10-'' (|Beers fc Christliebl [2005h . We 
must take care when comparing our results to observa- 
tions since we made the simplification that every Pop 
III star produces a PISN (for a detailed se mi-analytic 
model of m etal enrichm ent, see iTumlinsonI [2006) . As 
discussed in iHelmi et al.l . the galactic halo may be com- 
posed of remnants of galaxies that formed from high-cr 
density fluctuations, and dwarf galaxies originate from 
low-tJ peaks. In this scenario, the objects (or its rem- 
nants) simulated here would most likely reside in galac- 
tic halos at the present day. If we attempt to match this 
metallicity floor of 10^"* Z© in the galactic halo, this re- 
quires ~8Mq of metals produced for every Pop III star or 
roughly one in ten Pop III stars ending in a PISN. More 
likely, the current nearby dSphs are hosted by larger dark 
matter halos than we have been able to simulate to date. 
Hence too simple chemical evolution inspired extrapola- 
tions may be premature. 

5. SUMMARY 

Radiative feedback from Pop III stars play an impor- 
tant role in shaping the first galaxies. We studied the ef- 
fects of this feedback on the global nature of high-rcdshift 
dwarf galaxies, using a set of five cosmology AMR sim- 
ulations that accurately model radiative transfer with 



Abel, T., Anninos, P., Zhang, Y., & Norman, M. L. 1997, New 

Astronomy, 2, 181 
Abel, T., Bryan, G. L., & Norman, M. L. 2000, ApJ, 540, 39 
Abel, T., Bryan, G. L., & Norman, M. L. 2002, Science, 295, 93 
Abel, T., Croft, R. C, k. Hernquist, L. 2001, in The Shapes of 

Galaxies and Their Halos, eds. P. Natarajan, Yale, p. 119 
Abel, T. & Wandelt, B. D. 2002, MNRAS, 330, L53 
Abel, T., Wise, J. H., k. Bryan, G. L. 2007, ApJ, 659, L87 
Ahn, K. & Shapiro, P. R. 2007, MNRAS, 375, 881 
Anninos, P., Zhang, Y., Abel, T., & Norman, M. L. 1997, New 

Astronomy, 2, 209 
Barkana, R., & Loeb, A. 2001, Phys. Rep., 349, 125 
Barkat, Z, Rakavy, G., & Sack, N. 1967, Phys. Rev. Lett., 18, 379 
Barnes, J., & Efstathiou, G. 1987, ApJ, 319, 575 
Becker, R. H. et al. 2001, AJ, 122, 2850 
Beers, T. C., & Christlieb, N. 2005, ARA&A, 43, 531 
Bertschinser. E. 1995, ArXiv Astrophysics e-prints, 

arXiv:astro-ph/9506070 
Bertschmger, Ji. 2u0i, ApJS, 137, 1 

Blumenthal, G. R., Faber, S. M., Flores, R., & Primack, J. R. 1986, 
ApJ, 301, 27 

Bond, J. R., Arnett, W. D., & Carr, B. J. 1984, ApJ, 280, 825 
Bouwens, R. J., k. lUingworth, G. D. 2006, Nature, 443, 189 
Bromm, V., Coppi, P. S., & Larson, R. B. 2002, ApJ, 564, 23 



adaptive ray tracing. Additionally, we focused on the 
metal enrichment of the star forming halos and their as- 
sociated star formation histories. Our key findings in this 
paper are 

1. Dynamical feedback from Pop III stars expel nearly 
all of the baryons from low-mass host halos. The baryon 
fractions in star forming halos never fully recover even 
when it reaches a virial temperature of 10'' K. The baryon 
fraction is reduced as low as ~0.05 with SNe feedback, 
three times lower than the cases without stellar feedback. 

2. Baryons on average gain angular momentum as they 
are expelled in feedback driven outflows. When it is rein- 
corporated into the halo, it increased the spin parameter 
up to a factor of 4 with radiative and SNe feedback. 

3. The accurate treatment of radiative transfer pro- 
duces a complex, multi-phase ISM that has densities and 
temperatures that can span up to 6 orders of magnitude 
at a given radius. 

4. Pair-instability SN preferentially enrich the IGM to 
a metallicity an order of magnitude higher than the over- 
dense filaments adjacent to the sites of star formation. 

5. Once a SN explosion cannot expel the gas in its host 
halo, the mean metalhcity fluctuates around IQ-^-^/so 
as there may be a balance between SN outflows, cold 
inflows, and contained SNe ejecta. 

We conclude that Pop III stellar feedback plays an in- 
tegral part to early galaxy formation as it determines the 
characteristics of the first galaxies. 
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